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We investigate the phase behavior of symmetric AB diblock copolymers confined into a thin 
film. The film boundaries are parallel, impenetrable and attract the A component of the diblock 
copolymer. Using a self-consistent field technique [M.W. Matsen, J.Chem.Phys. 106, 7781 (1997)], 
we study the ordered phases as a function of incompatibility \ an d film thickness in the framework 
of the Gaussian chain model. For large film thickness and small incompatibility, we find first order 
transitions between phases with different number of lamellae which are parallel oriented to the film 
boundaries. At high incompatibility or small film thickness, transitions between parallel oriented 
and perpendicular oriented lamellae occur. We compare the self-consistent field calculations to 
Monte Carlo simulations of the bond fluctuation model for chain length N = 32. In the simulations 
we quench several systems from \N — to \N = 30 and monitor the morphology into which 
the diblock copolymers assemble. Three film thicknesses are investigated, corresponding to parallel 
oriented lamellae with 2 and 4 interfaces and a perpendicular oriented morphology. Good agreement 
between self-consistent field calculations and Monte Carlo simulations is found. 
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I. INTRODUCTION. 



Amphophilic polymers are model systems for investigating mechanisms of self-assembly. Joining chemically distinct 
polymers - A and B - at their ends to form an AB diblock copolymer prevents macrophase separation of the two 



species. In order to reduce the number of energetically unfavorable interactions between distinct blocks in a melt, 
the molecules self-assemble into complex morphologies. The morphology is selected via a delicate balance between 
the free energy cost of the internal interfaces and the conformational entropy loss as the molecules stretch to fill 
space at constant density. The phase diagram in the bulk has been investigated in much detail as a function of the 
relative length of the blocks / and the incompatibility %iV 0-§|- The morphologies found in copolymer melts and 
copolymer/homopolymer mixtures M resemble the spatially structured phases of other amphiphilic systems (e.g., 

t-H , lipid/water mixtures). 

From a theoretical point of view, polymeric systems are particularly convenient for investigating mechanisms of 
self-assembly. Only a small number of parameters describe the system, i.e., the fraction / of A monomers in the 
diblock, the molecule's end-to-end distance R e and the incompatibility ^JV, where x denotes the repulsion between 
monomers of different species and N the number of monomers per molecule. In general, polymeric systems are well 
describable by self-consistent field theories using the Gaussian chain model [Q-Q . For a wide range of temperature the 
theory accurately calculates the excess quantities of the internal interfaces (e.g., the interfacial tension, the bending 
moduli, or the enrichment of solvent). The understanding of these interfacial properties makes polymers suitable 

*^J ' microscopic model systems for investigating the statistical mechanics of interacting interfaces. 

If the copolymers are confined into a thin film, the interactions with the boundaries influence the morphology and 
its orientation. Controlling the orientation of the spatial structure over rather large length scales is important for 
many practical applications [pL Consequentially, much experimental effort has been directed towards a tailoring of 
surface properties and studying the morphologies in confined geometry. Lambooy et al. M studied symmetric dPS- 
PMMA diblock copolymers between silicon wavers and found transitions between parallel oriented lamellar phases 
with different numbers of internal interfaces upon increasing the film thickness. Kellogg and co-workers |Hj] studied 
a similar system between surfaces, which were coated with random copolymers. This boundary mimics neutral walls 
and for film thickness corresponding to 2.5 lamellar spacings in the bulk, indeed, perpendicular oriented lamellae 
could be observed. The perpendicular oriented morphology has attracted abiding interest as a template for lateral 
structure on the nanometer scale. 



The self-assembly was analyzed for strong segregation by Walton et al. 1 1 and in the framework of the self-consistent 
field theory by Pickett and Balazs |L2J and by Matsen fl3|,|l4| . Strictly neutral walls give rise to a perpendicular orien- 
tation of the lamellae. However, surfaces that interact favorably with one component stabilize parallel morphologies 
if the film thickness is compatible with the lamellar spacing in the confined state. Increasing the film thickness the 
self-consistent field studies revealed a sequence of perpendicular and parallel oriented lamellae. 



The orientation of the morphologies upon confinement makes thin films a promising candidate for investigating 
the details of their structure via Monte Carlo simulations. Copolymers in confined geometry have been first studied 
in Monte Carlo simulations by Kikuchi and Binder |15| . They found pronounced effects of the confinement on the 
ordering. However, similar to experiments, Monte Carlo simulations are plagued by very long relaxation times and 
reaching thermal equilibrium is rather difficult. Indeed, the mixed state of perpendicular and parallel orientated 
lamellae observed in the simulations and experiments |16[| was found to be unstable in framework of self-consistent 
field calculations |13j . Reviews of both experiments and theory on the self-assembly of block copolymers in thin films 
|L4|,[L7| give more details on previous work. 

The aim of the present work is twofold: On the one hand, we calculate the phase diagram of diblock copolymers in 
a thin film as a function of the film thickness and incompatibility. We use a self-consistent field technique developed 
by Matsen |0] . Both confining walls attract the A component of the diblock and we restrict ourselves to symmetric 
diblock copolymers / = 1/2. The temperature spans the weak and intermediate segregation limit. On the other hand, 
we compare the stability of different phases in the self-consistent field calculations and in the Monte Carlo simulations 
of the bond fluctuation model |18) for chain length N = 32 and at incompatibility \N — 30. 

Our paper is arranged as follows: First we briefly describe the self-consistent field technique for diblock copolymers 
in thin films and introduce the model used in the Monte Carlo simulations. Then we discuss the regions of stability 
of the various oriented lamellar phases as a function of the incompatibility and the film thickness. We present a 
qualitative comparison between the self-consistent field calculations and the Monte Carlo simulations and close with 
a brief discussion of our findings. 

II. MODEL AND TECHNIQUES 

We consider n diblock copolymers containing N — Na + Nb segments in a volume Ao x L x L. Ao denotes the 
film thickness, while L is the lateral extension of the film. The monomer number density in the middle of the film is 
denoted by p. The density at the film surfaces deviates from the density in the middle and it is useful to introduce 
the thickness A of an equivalent film with constant monomer density A = nN/pL 2 . The individual blocks of the 
diblock have identical length Na = Nb = fN, such that the diblock assembles into a lamellar phase in the bulk. R e 
denotes the end-to-end distance of the molecule. There is a short range repulsion between the two monomer species 
which can be parameterized by the Flory-Huggins parameter x- The two surfaces of the film are impenetrable and 
hard. Therefore, there is no formation of islands or hole defects at the surface. However, the free energy of a confined 
film is also an important ingredient in understanding the static and dynamics of pattern formation in thin films with 
free boundaries Jl9f . Both walls attract the A component of the diblock and repel the B component via a short range 
potential. 

A. Self-consistent field calculations (SCF) 

The computational technique in our self-consistent field (SCF) calculations is very similar to the work of Matsen 
p3[. In a boundary region of width A w the total monomer density drops to zero at both walls. In our calculations 
we assume the monomer density profile p$o to take the form 

< x < A w 

A w < x < A - A w (1) 

A - A w < x < A 

The width and the shape of the density profile near the wall is determined by a competition of the entropy loss of 
the polymers near the wall (favoring a thick boundary) and equation of state effects, which try to restore a spatially 
homogeneous density. The particular choice of the density profile is employed for computational convenience. The 
value of A w and the shape of the profile has little effect on the relative stability of the different phases. We choose 
A w = 0.15-Re in accord with the previous study [p~3f . This value is close to the ratio of the interaction range and the 
end-to-end distance in the Monte Carlo simulations. A film with the same number of monomers but uniform density 
would have the thickness A = Ao — A w . 

Both walls attract the A component of the diblock and repel the B component via a short range potential. The 
monomer wall interaction H is modeled as [Eo|: 




^^{l + cos(^)} for 0<a:<A w 

H(x) = { for A w < x < A - A w (2) 
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[l + cos ( " ( l° i ; x) ) } for A - A„ < a; < A 



The normalization of the surface fields Ai and A2, which act on the monomers close to the left and the right wall, is 
chosen such that the integrated interaction energy between the wall and the monomers is independent of the width 
of the boundary region A„, . 

The microscopic monomer densities <fr A and <£> B can be expressed as a functional of the polymer conformations 
K(r)}: 

M*) = -r£ / dr5(r-r Q (r)) (3) 

where the sum runs over all n diblock copolymers in the system and < r < 1 parameterizes the contour of the 
Gaussian polymer. A similar expression holds for <fr B (r). With this definition the partition function of a melt of 
Gaussian diblock copolymers takes the form: 

Z ~JT>[ T ]V[ T ] exp (-pj d 3 r { x *a*b - H(r)($ A (r) - $ B (r))}) 

x 6 f* (r) - *x(r) - &b(t)) (4) 



The functional integral T> sums over all chain conformations of the diblock copolymers and V[r] ~ 
exp I ~ 2 Nb i lo ^ T (jr) ) denotes the statistical weight of a non interacting Gaussian polymer, b 2 = R 2 /{N — 1) 
is the statistical segment length of the polymer. This simple model neglects the coupling between the interaction 
energy and the chain conformations |2l| ]. Thus, it cannot reproduce the stretching of the diblock copolymer in the 
disordered phase observed in simulations |22] and experiments p3|-p5||. Moreover, the Gaussian chain model neglects 
a finite stiffness of the polymers, and the chain extensions parallel to the walls (in the parallel orientated lamellae 
or disordered state) always remain unperturbed. The Boltzmann factor in the partition function incorporates the 
thermal repulsion between unlike monomers and the interactions between the monomers and the walls. The last 
factor represents the incompressibility of the melt in the center of the film and enforces the monomer density to decay 
according to Eq.(pl) at the walls. A finite compressibility of the polymeric fluid is neglected. 

Introducing auxiliary fields Wa, Wb, &a, &b and 5 we rewrite the partition function of the multi-chain system in 
terms of the partition function of a single chain 



•nil- nn- -r,«i\ rvi, -n- / ^\WA) Wb, $A> <Ab) ^j \ 

T>WaT>W b T>®aT>®bT > ^ exp ( — — ) (5) 



The free energy functional has the form: 

F[W A ,W B ,$ A ,<f> B ,E] _ 

T^, = - hi Q[W A , W B \ 

+ - f d 3 r X N$ A {r)$ B (r) 

- 1 J d 3 r H(r)N {$ A (r) - $ B (r)} 

- i / d 3 r {W A (r)<P A (r) + W B (r)$ B (r)} 

- 1 J d 3 r 3(r) {$ (r) - $ A (r) - $ s (r)} (6) 
where Q denotes the single chain partition in the external fields W A and W B : 

Q[W A , W B ] = ^J V[v\P[v] exp I - J dr W a (t(t)) - J At W e 



Mr)) (7) 



The functional integration in Eq.(|5|) cannot be carried out explicitly. Therefore we employ a saddlepoint approx- 
imation, which replaces the integral by the largest value of the integrand. This maximum occurs at values of the 
fields and densities determined by extremizing T with respect of each of its arguments. These values are denoted by 
lower-case letters and satisfy the self-consistent set of equations: 

w A (r)= X Ncj)B-H(r)N + ^r) 
w B (r)=xN<t>A + H(r)N + £{r) 
rf 



4>a{v)^Jw J At <5(r-r(r)) 

exp I - / At w a (y(t)) 
4>b{*) = t; f W J At S(t-t(t)) 
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At this stage fluctuations around the most probable configuration are ignored. Hence, the interfaces in the self- 
consistent (SCF) field calculations are ideally flat and there is no broadening by fluctuations of the local position of 
the interfaces (capillary waves). 

The distribution of copolymer segments is calculated by solving appropriate diffusion equations for the end segment 
distributions. In order to study the self-assembly into various morphologies, we expand the spatial dependence of the 
densities and fields in a set of orthonormal functions that possess the symmetry of the morphology being considered. 
This results in a set of non-linear equations which are solved by a Newton-Raphson like method. Substituting the 
saddlepoint values of the densities and fields into the free energy functional Q6J) we calculate the free energy of the 
various morphologies. For the perpendicular oriented lamellar phase the value of the free energy has to be minimized 
with respect to the lamellar spacing. We use up to 220 basis functions for the calculation of the perpendicular oriented 
lamellar phases and thus achieve a relative accuracy of the order 10 -4 for the free energy. 

B. Monte Carlo simulations (MC) 

For the Monte Carlo (MC) simulations we employ the bond fluctuation model Pq|. This coarse grained lattice model 
captures the relevant universal features of polymeric materials: excluded volume of segments, chain connectivity, and 
a short range thermal interaction. Many thermodynamic properties of the model have been determined in previous 
studies [£6| and the model is a good compromise between the computational advantages of a lattice model and a 
faithful representation of continuum space properties. In particular, the relation between the model parameters, the 
local fluid-like packing of monomers, and the phase behavior has been investigated in detail. Within the framework of 
this model a small number of chemical repeat units is represented by the eight corners of a cube on a three dimensional 
lattice. Monomers along a polymer are connected via one of 108 bond vectors of length 2, \/E, \/6, 3, and yTO. The 
distances are measured in units of the lattice spacing. The polymers comprise N — 32 monomers. We work at 
a monomer number density p = 1/16. This value corresponds to a concentrated solution or a melt. Under these 
conditions the end-to-end distance is R e « 17 and, in accord with previous studies, we use b = R e /\/N — 1 = 3.05 for 
the statistical segment length in the SCF calculations p7[ . 

One half of the polymer consists of A monomers, the other consists of B monomers. Monomers of the same type 
attract each other via a square well potential, while there is a repulsion between unlike species. The interaction 
range comprises the 54 nearest neighbors up to a distance Vo. This corresponds roughly to the first neighbor shell 
in the monomer density pair correlation function. The well depths of the interactions are chosen symmetrically: 

£AA — (BB = —CAB = -(■ 



The phase behavior and structure of a binary blend of A and B homopolymers 1 28 29] has been investigated in the 
framework of this model. The phase diagram of binary homopolymcr mixtures and ternary homopolymer/copolymer 
blends |30| as well as the interfacial structure between unmixed phases are well describable by the Gaussian chain 
model if the Flory-Huggins parameter is identified via |3lJ ] 

h '" d 3 r 5 intcr (r) « -^- 2.65 , (9) 



ksT J r <^/6 ksT 



where we have used the symmetry of the monomeric interactions. Here, g lnteT denotes the intermolecular pair correla- 
tion function in the melt and the integral is extended over the range of the square well potential. Moreover, we have 
assumed that g lntcr is largely independent from temperature. We use the value 2.65 for the number of monomers of 
different chains in the range of the square well interaction. This identification of the Flory-Huggins parameter is based 
on the energy of mixing; entropic contributions to the free energy due to packing effects or conformational changes 
are negligible. We perform Monte Carlo (MC) simulations at e = 0.1769fcsT. This value corresponds to x^ = 30. 

The simulation cell possesses a geometry of the form Ao x L x L. Periodic boundary conditions are applied in 
the lateral directions. The two impenetrable walls at x = and x = Ao + 1 interact with monomers in the nearest 
two layers via a square well potential. An A monomer in the interaction range of the walls lowers the energy by an 
amount e w while a B monomer increases the energy by the same amount. If the surface is completely covered with 
A monomers the wall interaction energy per chain is: 

F wall = _ 2 2e^V 

ALT ^ ; 



nksT A/c^T 

In the simulations we use the value e w — O.lfcsT. Simulations p2| of binary blends (with identical interactions) show, 
that for this strength of surface interactions the A component wets the surface for e < 0.043fcsT or %iV < 7.3. For the 
comparison between the SCF calculations and the MC simulations we adjust the surface fields A as to result in the same 
contribution to the energy. Hence, the surface interactions in the MC simulations correspond to A\N = A2 N = 0.375 
(cf. Eq.pj) below). 

The MC simulations comprise three different moves: The conformations of the polymers are updated via local 
hopping of the monomers and slithering snake like motions. In the former, we randomly choose a monomer and 
try to displace it by one lattice unit in a random direction. During the slithering snake attempts, we randomly 
choose a chain end and try to attach it at the opposite end of the chain. The monomer identity (A or B) of the 
chains is correspondingly updated to precisely conserve the composition of the chain p2| . The latter moves relax the 
conformations of the polymers a factor of the order N faster than local updates. Moreover, we allow for A ^ B flips, 
in which the identity of the A and B monomers of a randomly chosen polymer are exchanged. One Monte Carlo step 
consists of 3 slithering snake attempts per chain, 1 local hopping attempt per monomer and 1 A ^ B flip per diblock. 
Every 12500 Monte Carlo steps a configuration was stored for further analysis. Since we are interested in studying 
the stability of different morphologies, we do not impose a specific morphology on the starting configurations [J33[ . 
Rather, we let the structure self-assemble via a quench from the disordered phase to the ordered state at xN = 30, 
and monitor the morphologies which occur in several independent runs with identical parameters. 

III. PHASE DIAGRAM. 

At high incompatibility, many features of confined diblock copolymers can be deduced from the strong stretching 
theory M. This has been applied to study the effect of confinement by Turner fl34] and Walton et al fO] (see also 
fL5|). We follow the notation of Ref. fl3[| . We consider a parallel lamellar phase L p with p internal interfaces. In 
the limit \N 3> 10 the A and B rich domains are well segregated and the junction points are confined to a narrow 
interfacial region. To fill space uniformly the copolymers stretch. In a lamellar morphology each half of the diblock 
forms a brush. These brushes do not interpenetrate. In the parallel lamellar phase L p with p interfaces each brush 
has the height A/2p and the free energy cost due to the stretching of the chains in the brush |3^] amounts to: 

^brush _ ^ 2 (A/2 P ) 2 

nk B T 8{N/2)b 2 [ ' 

Each half of the diblock covers an interfacial area pN/pA. Estimating the value of the interfacial tension between A 
and B domains by the interfacial tension in a binary blend a/ksT « p&\/x/6, the free energy contribution of the 
internal interfaces per polymer takes the form: 

r ' '•'•" (12) 



nksT 

The balance between these two terms determines the behavior in the bulk. This leads to a prefered lamellar spacing 
£>b = 2A m ; n /p = 2(8xA/37r 4 ) 1 / 6 i? e . If the walls preferentially interact with one component the interaction energy 
of the monomers with the walls gives another contribution to the free energy. Using the expression for the density 
profiles and the wall monomer interactions, this contribution takes the form: 



wall 



, T 4/d'rg(r) ( tf,(,) s - A ' y y (13) 

where we have assumed that each surface is covered completely with the energetically favored component. For parallel 
lamellar phases with an odd or an even number of interfaces exposed to symmetric or antisymmetric surfaces fields, 
respectively, the contributions cancel. 

The confinement into a film also reduced the conformational entropy of the molecules. Using ground state domi- 
nance, the entropy of an inhomogeneous melt is given by pq | : 

^o„f _ 6 2 /■ (V$ (r)) 2 ttW 

A („\ O/l A A V / 



nk B T 2An J $ (r) 24A W A 

This contributions depends strongly on the detailed density profile at the wall. However, it does not discriminate 
between the different phases and, hence, is irrelevant to the stability of the different morphologies. 

In the perpendicular lamellar phase Lj_ both components cover an equal amount of surface area and, hence, 
the contribution of the surface fields vanishes. Moreover, the lamellar period is free to adjust as to minimize the 
contribution from the internal interfaces and the chain stretching. As a result, the lamellar spacing in the strong 
stretching theory is identical in the perpendicular morphology of a film and in the bulk. 

FigjI] presents the excess free energy for the different phases in the strong stretching approximation for \N = 30. 
We plot the difference (F — F\,)A/nkBTR e , which is proportional to the difference between the free energy in the 
confined geometry and the bulk free energy per unit area of the film. For neutral walls (a), which do not prefer a 
component, the L± phase is stable for all film thicknesses. The parallel lamellar morphology has the identical free 
energy only if the film thickness coincides with half integer multiples of the bulk period D\>. For symmetric walls 
(b) the free energy of the parallel morphology with an even number of interfaces is lowered, because the prefered 
component is brought into contact with both walls. The free energies of the other phases remain unchanged. Upon 
increasing the film thickness one observes transitions from morphologies with an even number 2p of interfaces at film 
thickness around pD^ and perpendicular orientated lamellae for thickness (2p + l)D\,/2. In the case of antisymmetric 
surface fields (c) the free energy of the parallel lamellae with an odd number (2p+l) of interfaces is lowered, and 
one finds transitions between those parallel lamellae for film thickness close to (2p + l)£>b/2 and the perpendicular 
morphology for thickness around integer multiples of the bulk period. 

In Fig we compare the results with the full SCF calculations. In this representation the curves are qualitatively 
similar. In particular, both approaches yield the same sequence of morphologies as the film thickness is increased. 
However, the absolute values of the free energies differ by more than 20% for these parameters. Moreover, there are 
some subtle differences: For neutral walls (a), the free energy of the perpendicular morphology is strictly lower than 
that of the parallel oriented lamellae. This is in accord with the calculations of Pickett and Balazs JL2]. For symmetric 
surface fields (b) the free energy of the parallel morphology with an even number of interfaces is lowered. The value 
of the relative shifts upon increasing the surface field in the strong stretching theory and the SCF calculations agree 
nicely. In both approaches the free energy of the parallel morphology with an odd number of lamellae remains almost 
unaffected. This indicates that the structure at the surfaces is hardly perturbed by the weak surface fields. However, 
the free energy of the perpendicular morphology is independent of the surface fields in the strong stretching theory, 
while the presence of the surface fields lowers the free energy in the SCF calculations, indicating a dependence of the 
spatial arrangement on the surface fields. 

The composition profiles of the perpendicular morphology are presented in Fig.H for neutral walls (a) and symmetric 
walls attracting the A component (b). The A rich regions are bright and B rich regions are shaded darkly. In the 
case of neutral walls, the interface between the A and B domains runs strictly perpendicular to the surface. However, 
the interfacial width broadens close to the surface. Partially, this is due to the reduction of the density in the surface 
region, which reduces the incompatibility between the two components. This reduces the AB interfacial tension as 
the interface intersects the wall. Moreover, the polymers in the vicinity of the surfaces are aligned parallel to the wall, 
and this orientation is compatible with their conformation at the AB interface as it approaches the wall JL2J. Both 
effects reduce the free energy of the AB interface in the perpendicular morphology. This gives rise to a negative line 
tension |13|] and tends to stabilize the perpendicular phase. 

Upon increasing the surface interactions (b), the AB interface bends and intersects the wall at an angle. This 
distortion of the interface close to the surface increases the surface area covered by the energetically favored A 
component and lowers the free energy of the perpendicular phase |12fl. This effect is not captured by the strong 
stretching approximation. However, Pereira and Williams [B7J have argued that it remains typically small. Note that 
the surface fields are rather small such that the A component does not wet the surface. 

The phase diagram as a function of the incompatibility \N and film thickness is presented in Fig.0 for symmetric 
boundary fields. At high incompatibilities we find a sequence of perpendicular aligned lamellae L± and parallel 



lamellae L p with an even number of interfaces. The latter are stable for film thicknesses close to integer multiples 
of the bulk period. Upon increasing the film thickness the stability region of the perpendicular oriented morphology 
decreases. The free energy difference between the two morphologies is related to the balance between the surface 
interactions which favor parallel orientation and the free energy costs of imposing a lamellar spacing which differs 
from the prefered bulk value. The surface contribution is independent of the film thickness. The free energy costs 
per lamellae due to a mismatch in the lamellar spacing increase quadratically. Hence, it is favorable to distribute the 
mismatch evenly among the p lamellae; the mismatch per lamellae is proportional to 1/p. Therefore the free energy 
of the film due to deviations of the film thickness from the prefered spacing decreases like p x (1/p) 2 ; i.e., for thick 
films the mismatch becomes unimportant and only parallel lamellae are stable. 

In thin films the translational symmetry parallel to the walls is spontaneously broken and the SCF theory predicts 
a second order transition from the disordered state to the perpendicular lamellar phase. For thicker films we find 
direct transitions between the parallel lamellar phase with 4 and 6 interfaces. Only at higher incompatibilities we 
encounter a triple point, at which the two parallel lamellar phases L4 and Lq coexist with a perpendicular aligned 
lamellar phase L±. At higher incompatibilities we find the sequence L^Lx and Lg upon increasing the film thickness. 
We expect this behavior to be representative for larger film thicknesses and the incompatibility at which the triple 
point occurs increases with the film thickness. An important point is that the theory predicts a gradual onset of 
parallel ordering (L p ) as \N increases, without a phase transition from the disordered phase. This happens because 
for finite A the surface fields create surface induced order of lamellar type already in the disordered phase, and this 
order gets gradually stronger as \N increases. 

The composition profiles in a thin film close to these critical points, where the perpendicular lamellae emerge, 
are presented in Fig|| for a film thickness A/R e = 0.55 and Fig.|| for A/R e = 1.92. When the incompatibility is 
increased the perpendicular modulations become more pronounced. Slightly above the critical points the i?-rich 
domains for A/R e = 0.55 and the A-rich domains for A./R e = 1.92 form cylinders which run parallel to the surfaces. 
This behavior resembles the fingerprint-like morphology observed in experiments of Chaikin and co-workers [ p8[ . For 
slightly asymmetric diblocks or stronger surface fields even more pronounced effects could be anticipated. However, 
the neglect of fluctuations imparts a quantitative inaccuracy to the SCF calculations in the weak segregation limit. 
In particular, the existence of critical points where a second order transition from the disordered phase into the 
perpendicular oriented lamellar structure occurs is questionable. In the bulk, one encounters a fluctuation-induced 
first order transition rather than a critical point [ p9| . 

IV. COMPARISON BETWEEN SELF-CONSISTENT FIELD (SCF) CALCULATIONS AND MONTE 

CARLO (MC) SIMULATIONS. 

In the following we compare our SCF calculations at \N = 30 and K\N = A2A = 0.375 to the corresponding MC 
simulations in the framework of the bond fluctuation model. The value of the incompatibility lies in the regime of 
intermediate segregation and is well inside the experimentally accessible range. For much smaller incompatibilities 
there are strong fluctuation effects. We have performed some preliminary simulations in a 256 x 64 x 64 geometry 
with e w — O.lfcsT to study the ordering behavior. The results for the difference between the A and B monomer 
density in the vicinity of the surfaces are displayed in Fig.0. Upon increasing the incompatibility, the amplitude 
and correlation length of composition fluctuations increase. Moreover, the period of the oscillations increases slightly 
upon increasing the incompatibility. This indicates a stretching of the molecules. For e = 0.09 we observe a weak 
modulation of the composition across the whole film thickness. Hence, we expect the order-disorder transition to 
occur in the range e « 0.09(1) or 13.5 < %A < 17 J4(|. This is in accord with a previous estimate of the transition 



temperature |30|. Moreover, the magnitude of this shift in the transition temperature is in qualitative agreement 
with fluctuation corrections calculated by Fredrickson and Helfand p9[ . They predict XcN = 10.5 + 41/(_R 2 /r 2 / 3 ) 1 / 3 
for very long chains. Deviations of similar magnitude, though by a different mechanism, have been predicted in the 
framework of the P-RISM theory []41[ . For chain length N = 32 we anticipate rather large deviations from the SCF 
calculations at incompatibilities smaller than xA = 25. Since for the present model the transition point is not known 
to high precision, we did not attempt a quantitative co mp arison of the results in Figj^ with theoretical results for the 
order parameter profiles for surface-induced ordering fl42fl . 

At stronger incompatibilities ee > 0.5, there occurs a phase separation into a homopolymer-rich phase and a phase 
rich in vacancies. This would correspond to %A ss 85 according to Eq.(|9|). However, we do not expect this equation 
to hold in this temperature regime, because the fluid structure at these temperatures differs significantly from the 
high temperature structure. At much smaller incompatibilities, already, the width of the internal interfaces becomes 
comparable to the length scale of the local polymer architecture. Previous studies have shown that in the strongly 
segregated regime the Gaussian chain model may yield qualitatively erroneous results p8J. These two considerations 



yield a rough estimate for the temperature interval in which good agreement between the MC simulations and the 
SCF calculations can be expected. 

We have quenched three different systems from their athermal state (e = 0) to e = 0.1769/csT. The systems have 
the geometry 30 x 96 x 96, 46 x 93 x 93 and 56 x 96 x 96 in units of the lattice spacing. According to the SCF calculations 
these geometries correspond to the L 2 , the Lj_, and the L4 phase, respectively. The lateral extension in the Lj_ phase 
has been chosen compatible with the lamellar spacing in the SCF calculations. For each of these geometries we have 
simulated at least 4 independent systems. We do not find transitions from one morphology to a different one during 
the simulation time which exceeds at least 3.8 • 10 6 Monte Carlo steps for each system. Hence, we cannot rule out 
metastability effects completely. However, the simulation time is long enough that the observed structures are free 
of defects (on the scale of the simulation cell) and the composition profiles of systems which assembled into the same 
morphology agree (cf. also El). 

In Fig.g we present the monomer density profiles of the L 2 phase in the SCF calculations (a) and in the MC 
simulation (b). Qualitatively, the profiles are similar. Both data sets show almost completely segregated A and B 
rich regions separated via two interfaces. However, in the SCF calculations the total density rises smoothly from zero 
at the film boundaries to one in the middle of the film according to Eq.(|l|). In the MC simulations the monomers 
pack against the wall and produce oscillations in the density profile near the surfaces. These details of the local fluid 
structure are not captured in the Gaussian chain model. Note that the average density in the two layers nearest to 



the walls is close to the bulk density. Thus the simple estimate for the energy contribution of the walls (cf. Eq.(10)) 
is largely unaffected by the packing effects. 

Snapshots of the final morphologies are presented in Fig.|[ In panel (a) we present the final snapshots for film 
thickness A = 30. The surfaces correspond to the top and bottom plane. All systems have assembled into the L 2 
phase. In the snapshots the A component corresponds to the darker species. Half an A lamella is located at each 
wall. However, one also observes rather strong fluctuations which allow the B component to protrude up to the wall. 
All films of thickness Ao = 46 (cf. Fig.pl(b)) have assembled into the L± phase. However, the simulations exhibit two 
different repeat distances. In the two systems, where the lamellae are oriented parallel to the box axis, the lamellae 
are spaced at a distance D = 1.827i? e , which agrees with the SCF calculations -Dscf = 1.822i? e - However, in the 
six other systems, the lamellae make an angle with the box axis and the repeat distance is D — 1.9387? e . This 
value exceeds the SCF prediction by 6%. In the SCF calculations such a deviation from the prefered repeat distance 
increases the free energy by 0.01 ksT per molecule or 8 fcgT for the whole system. Hence, the occurence of the 
larger spacing cannot be explained by thermal fluctuations alone. The difference might be traced back to the fact 
that we use for the SCF calculations the chain extension R e corresponding to the athermal state. As it has been 
observed in other simulations |2^] (cf. also Fig.Q) , the chains stretch even in the disordered phase to avoid energetically 
unfavorable contacts between the different blocks. For the comparison with the SCF calculations J44J we employ the 
two configurations in which the lamellae are oriented parallel to the box axis. Three films of thickness Ao = 56 (cf. 
0(c)) assembled into the £4 phase, whereas one system prefered the L± phase. For the comparison with the SCF 
calculations |44|] the latter system was discarded. Though we cannot rule out that some of the systems are trapped 
in metastable conformations it is very gratifying that we observe the morphologies predicted by the SCF calculations 
except for one case at the largest film thickness. 

V. SUMMARY. 

We have presented SCF calculations and MC simulations for symmetric diblock copolymers confined into a thin 
film. Both surfaces attract the same component of the diblock via a short range potential. We have calculated 
the phase diagram as a function of the incompatibility %iV and the film thickness in mean field approximation and 
discussed the stability of parallel L p and perpendicular L± aligned lamellar phases. At high incompatibility we find 
the sequence L 2 ,L^,L^L\^, L§ while we find a direct transition between the L4 and the Lq phase at weak segregation. 

At low incompatibilities we find rather pronounced deviations from the SCF theory. Most notably, for chain length 
N = 32 the onset of ordering occurs around 13.5 < XcN < 17 instead of Y^ CF iV « 10.5. The order of magnitude of 
the shift is compatible with corrections to the mean field behavior |jjj[llj. At very high segregation, we expect the 
local structure of the model to be important. 

At \N = 30 we have compared the results of the SCF calculations with MC simulations in the framework of 
the BFM. We find qualitative agreement between the MC simulations and the SCF calculations. In particular, we 
observed the L 2 , L± and L4 phases as predicted by the SCF calculations. The main difference between the SCF 
calculations and the MC simulations is the structure at the surfaces. While the density decays smoothly to zero in the 
SCF calculations, there are pronounced packing effects in the MC simulations. In both schemes the incompatibility 
at the film surface is reduced. In the SCF calculations it stems from the reduced density at the surface, while in the 



MC simulations it is due to the finite extension of the interactions ("missing neighbor effect"). Both effects give rise 
to a negative line tension |L3|] as the interface approaches the surface and lead to a stabilization of the L± phase. The 
effect is, however, more pronounced in the MC simulations. Another difference between the MC simulations and the 
SCF calculations is the dependence of the chain extension on the incompatibility. The majority of systems in the L±_ 
phase has assembled into a morphology with a lamellar spacing which exceed the prediction of the SCF calculations 
by about 6%. This goes along with a stretching of the diblock copolymers already above the ordering transition. 
The effect has been observed in previous simulations p2,E3| and experiments [p3|-p5[, and can be rationalized via a 
coupling of the intramolecular energy and the chain conformations. 

In view of these effects the agreement with the SCF calculations is satisfactory. A detailed comparison of individual 
profiles between the SCF calculations and the MC simulations shall be presented in the following paper jTJJ . 
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FIG. 1. Free energy as predicted by the strong stretching theory as a function of the film thickness at \N = 30. 
The solid lines correspond to parallel lamellar phases L p , where p denotes the number of AB interfaces. The 
dashed line represents the perpendicular lamellar phase, (a) neutral walls K\N — A2-/V — (b) symmetric walls 
AiiV = A 2 iV = 0.2, and (c) antisymmetric walls AiN = -A 2 N - 0.2 
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FIG. 2. Comparison between the strong stretching theory (dotted lines) and the self-consistent field theory 
(thick solid lines (L p ) and dashed dotted line (L±)) at \N = 30. (a) neutral walls AiiV = A2N — and (b) 
symmetric walls AiN = A2N = 0.2 
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FIG. 3. Contour plots of the composition in the perpendicular lamellar phase at \N — 30 and A/R e = 2.3. 
A rich regions are light, B rich regions are dark, (a) neutral walls AiiV = A2-/V = and (b) symmetric walls 
AiiV = A2A^ = 0.2 Lines correspond to constant compositions in the xz plane. The coordinate x across the film 
is oriented along the ordinate here. (Lengths are measured in units of R e .) 
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FIG. 4. Phase diagram of a thin film with symmetric walls AiA r = K2N = 0.2 as a function of the incom- 
patibility \N and the film thickness A/R e . L2, L4, and Le denote parallel lamellar phases with 2,4, and 6 AB 
interfaces, whereas L± denotes the perpendicular lamellar phase. The dashed lines mark multiples of the bulk 
lamellar period. The square denotes the approximate location of the triple point. 
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FIG. 5. Contour plots of the composition in a symmetric film A/R e — 0.55 and AiiV = A2N = 0.2. Upon 
increasing the incompatibility, one encounters a second order transition to a perpendicular lamellar phase. The 
critical point occurs at about x^N — 11.72. (a) \N = 11.72, (Figure available upon request) (b) \N = 11.8, and 
(c) X N = 12.5 
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FIG. 6. Contour plots of the composition in a symmetric film A/R e — 1.92 and AiiV = A2N = 0.2. Upon 
increasing the incompatibility, one encounters a second order transition to a perpendicular lamellar phase. The 
critical point occurs at about XcN = 10.87. (a) x^ = 10.85, (b) x^ = H-5, and (c) \ N = 13 
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FIG. 7. Difference between the A and B monomer density in the vicinity of the surface. Upon approaching 
the order-disorder transition, we observe that the correlation length and amplitude of composition oscillation 
increase. 
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FIG. 8. A and B density profiles of the L2 phase in a symmetric film A/R e — 1.71 and AiiV = A2iV = 0.375. 
(a) results of the self-consistent field theory and (b) results of the Monte Carlo simulations. Note the packing 
effects at the wall. 



FIG. 9. Snapshot of the morphologies found in the MC simulations by quenching the systems from xN — 
to \N = 30. (a) Aq = 30 (b) Aq = 56 (c) Ao = 46 Figures available upon request. 
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